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Abstract. The Airy transform is an ideally suited tool to treat problem in 
classical and quantum optics. Even though the relevant mathematical aspects 
have been thoroughly investigated, the possibility it offers are wide and some 
aspects, as the link with special functions and polynomials, still contains un- 
explored aspects. In this note we will show that the so called Airy polynomials 
are essentially the third order Hermite polynomials. We will also prove that 
this identification opens the possibility of developing new conjectures on the 
properties of this family of polynomials. 



1. Introduction 

The theory of ordinary and generalized Hermite polynomials has largely ben- 
efited of the operational formalism. The two variable Hermite-Kampe de Feriet 
polynomials [1] can be defined using the following identity [2]: 



(1.1) H n (x,y) = 



x 



which involves the action of an exponential operator, containing a second order 
derivative, on a monomial. The explicit form of the polynomials H n (x,y) can be 
obtained by means of a straightforward expansion of the exponential in eq. ( 1.1 1, 
which yields: 



[n/2] 



„n—2r , 



d- 2 ) H ^y) = E fr % r * n = " ! E ' 

r=0 ' r=0 V '' 

where the variables x, y independent each otheiQ By keeping, therefore, the de- 
rivative of both sides of eq. (1.1 1 with respect to y, we find that the Hermite 
polynomials can be viewed as the solution of the following heat equation 

(1.3) d y F(x,y) = d 2 x F(x,y) , F(x,0) = x n . 

For y > they can be written in terms of the Gauss- Weierstrass transform [3] : 

(1.4) H n (x,y) = f°° d£r exp{-^-^ 

which is a standard mean of solutions for the heat type problems. 

The higher order Hermite polynomials [2], widely exploited in combinatorial 
quantum field theory can be expressed as a generalization of the operational 



By interpreting the variable y as a parameter, the standard Hermite form are recovered by 
the identities H n (2x, —1) = H n (x) and H n (x, —1/2) = He n (x). 
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identity (1.1 1, and indeed they writ^] 

[n I m] 

(1.5) HW(x,y) = e yd * x n = nl 



x n~mr y r 



— ' (n — mr)\ r! 

Therefore, we can ask whether an integral transform, a sort of generalization of 
the Gauss- Weierstrass transform, also holds for the higher order case. We start 
discussing the case of Hermite polynomials of even order with negative values of 
the y parameter, namely: 

(1.6) HW>(x,-\v\) = e-M*? x n , 

We express this family of polynomials in terms of a suitable transform following 
the procedure, put forward in [5], which considers the operator function 

(1.7) F = f(d x ) , 

where f(x) is a function admitting a Fourier transform f(k). With this assumption 
we find that the operator F can be written as: 

1 f°° 

(1.8) F = —f== j d*/(*)e** 8 « , 

and, therefore, we can express the action of the operator Fona given function g{x) 
as the integral transform indicated below 

(1.9) Fg(x) = — dkf(k)e* ka *g(x) = — / dfc f(k) g(x + tk) . 



We can now apply the same procedure to express the exponential operator inter- 
vening in the definition of iJ«, (see eq. ( 1.5 )), thus obtaining the following integral 
transform yielding the even order Hermite polynomials 

, , i r°° / , \ n 

(1.10) Hto\x,-\y\) = -= / dke 2p (k) (x-ik 2 V\y\) 



/2tt 

with 



e 2p (k) = I dxe x ~ 

v27T 



ikx 



whicrj^J after a redefinition of the variable, can also be written as 

1 1 f 00 
(1.11) H^(x,-\y\) = d^"e 2p 

V27T 1 2 {/\y\ J-oo 




W\ 

It must remarked that the same procedure, applied to the case m = 2, does not 
lead to a Gauss- Weierstrass transform as in eq. (1.4 1, which holds only for y > 0. 

In the next section we will extend the formalism developed in these introductory 
remarks, and prove that the Airy transform and the associated polynomials can be 
framed within the same context. 



The upper index m, denoting the order of the polynomials, is omitted for Hermite polynomials 
of order 2. 

■^The choice of discussing even order m is motivated by the the request that the integrals 
defining e m (fc) is convergent. 
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2. The Airy transform and the Airy polynomials 
The higher order Hermite polynomials satisfy the following recurrences [2J: 
d x Hl?\x,y) = nH^\(x,y) 

(2.1) H%\(x,y) = xHt\x,y)+m- ^—H^l +1 (x,y) 

[n — m + lj! 

d y H^(x,y) = ^^Hjri m (x,y). 

with the combination of the first and third recurrences yielding 

(2.2) d y H^ (x, y) = d?H<T> (x, y) , (x, 0) = x n , 

Therefore, the higher order Hermite polynomials satisfy a generalized heat equation, 
and this justifies the operational definition given in eq. (1.5). Furthermore, by 
interpreting y as a parameter, we can use the first two recurrences to prove tha 
they satisfy the m-th order ODE: 

(Am j \ 

y^ + x dx) H ™ {x > v) = nH " m){x > v) 

In the previous section we have considered even order Hermite polynomials only. 
Here we will discuss the third order case and their important relationship with the 
Airy transform and the Airy polynomials [6] . Before getting into this specific aspect 
of the problem, we note that the following identity holds: 

^Xx 3 / ii „s/\xt 



(2.4) e Xx = / dte VAxt Ai(t) , 

J —OO 

where 

(2.5) Ai(t) = i- J°° d £exp{^ 3 + ^} , 

is the Airy function [7]. 

Let us now consider the third order PDE 

(2.6) d y F(x,y) = dlF{x,y) , F(x,0) = f(x) 
whose formal solution writes 

(2.7) F(x,y) .■""•/(,•). 



By applying the identity given in eq. ( |2.4[ ), and by limiting ourselves to the case 
y > 0, we find: 



F(x,y) = ± jH dfcAz e^ kd * f(x) 

dkAi(J^=j f(x + yyk) 



J_ f°° ,, .. ( k 



OO 
OO 



which is recognized as the Airy transform of the function f(x). The concept of the 
Airy transform was introduced in ref. [7J , and has found noticeable applications in 
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classical and quantum physics [6] . The Airy transform of a monomial, namely 



(2.9) 



i n (x,y) = 



1 



/3y 



d£,Ai - 



/3y 



(y>o) 



has been defined as Airy polynomials [HE], but, according to the discussion of the 
previous section, they are also recognized as the third order Hermite polynomials 
Hn '(x, y). The characteristic recurrences (2.1 ) (specialized to the case m = 3), can 
also be directly inferred from eq. (2.9). 



3. Concluding remarks 

For further convenience we will introduce the following two variable extension of 
the Airy functions 

1 [°° 

(3.1) Ai(x, y) = — / d£ exp {iyt 3 + ixt) , 

2 71 " J -co 

that, in terms of the ordinary Airy function writes 

(3 - 2) Ai ^ y) = wv A iwv 

It is also easily shown that it satisfies the ODE 

d 2 

(3.3) 3y-^Ai(x, y) - xAi(x,y) = 

and that any other function linked to Ai{x,y) by 

(3.4) AiS m \x,y,z) = e 29 ™ Ai(x,y) 
satisfies the ODE 

(3.5) 3y^rAi^(x,y,z) - e zd * [xAi(x,y,z)} = 

ax z 

which on account of the identity 

(3.6) e zd * x = (x + mz^r 1 ) e zd ™ 
can be written as 

(3.7) \mz-^^-3y^ + xjAi^ m \x,y,z)=0 

If we assume m = 2p + 1, y = 0, and mz = 1, we get the following generalization 
of the Airy function 

(3.8) Ai( 2p+1 Hx) = — [ dt exp I t 2p+1 +ixt 

2tt J_ oc y 2p + 1 

satisfying the ODE 

(3.9) Q^ + ^^+i) (a0 = . 

As an example a comparison between the ordinary Airy function and its general- 
ization of order 7 is shown in Fig. [l] 

The point of view developed in sec. 2 can be generalized to the case with m > 3. 
We consider the generalized Hermite polynomials of odd order (m = 2p+l) and note 
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Figure 1. Comparison between the seventh order Airy function 
(red curve) and the ordinary Airy function (blue curve). 



that the same procedure exploited in the previous sections leads to the following 

integral representation 

(3.10) 



H^(x-\y\) = 



a "V(2p + i)M J- 



de^ (2p+1) 



^(2p+l)\y\ 



The use of a further variable or parameter in the theory of special functions 
and/or polynomials offers a further degree of freedom, which may be helpful to 
derive new properties, otherwise hidden by the loose of symmetry deriving from 
the fact that a specific choice of the variable has been done. This is indeed the case 
of the Hermite polynomials, which acquire a complete new flavor by the use of the 



y variable and the case of the Airy function given in eq. (3.1 ) which is easily shown 
to be the natural solution of the PDE 



(3.11) d y Ai(x,y) 
and the translation property 



-8%Ai{x, y) 



(3.12) 



-z9t 



Ai(x,y) = Ai[x, y + z) 



According to eq. (3.11), we can write the solution of the equation 
(3.13) d y F(x,y) = -d 3 x F(x,y) , 

as follows 



F(x,0) = g(x) , 



(3.14) 



/CC 
dZAi(x-Z,y)g(Q. 
-cc 
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It is evident that the further generalization 

(3.15) Ai^ p+1 \x,y) = — J dt exp{iyt 2p+1 + ixt) , (y > 0) , 

satisfies a higher order heat equation and can be exploited to obtain the solution 
of the same family of equation in terms of an appropriate integral transform. 
Before concluding this note let us consider the following Schrodingcr equation 

h 2 

(3.16) ih$(x,t) = -— dl^{x,t) + FxV(x,t) , *(z,0) = ip(x) , 

2m 

describing the motion of a particle in a linear potential (F is a constant with the 
dimension of a force). This equation can be written in the more convenient form 

(3.17) id T ^(x,r) = -dl^{x,r) + bx^{x,r) , *(ac,0) = ip(x) , 
wher^3 

ht , 2Fm 

As discussed in refs. [8], the previous equation can be treated using different 
means, a very simple solution being offered by the use of the following auxiliary 
function: 

(3.19) <f>{x,r) = e^ 9 * *(x,r) 
which satisfies the equation 

(3.20) i3 t $(i,t) = 6x$(x,r) . 



The solution of eq. (3.17) can therefore be written as 

(3.21) $(x,t) = Aifj{x) , A = e - ^e- ,6BT e^ , 

where the operator A, once written in an integral form, is recognized as the Airy 
transform. 

The solution of the equation 

(3.22) id T V{x,T) = -dl P V{x,T) + bx^(x,T) , *(x,0) = ip{x) , 
can be obtained in an analogous way, and it is given by 

(3.23) V(x,t) = A^ip(x) , i (2p) = e"^ 8 - ,+1 e-' kT e(W 8 ' f+1 , 

where the operator A^ 2p ^ can be expressed in terms of integral transforms, linked 
to the higher order Airy functions discussed in these sections. Further examples 
of generalization of the Airy function are due to Watson [5], who introduced a 
complete class of functions with interesting properties for physical applications. As 
an example we consider the function expressed as 



POO 

(3.24) W(x) = dt cos(i 4 + 4xt + 2x 2 ) 

Jo 

which is compared to the ordinary Airy function in Fig. [2j and satisfies the ODE 

d 2 

(3.25) W(x) + 4x 2 W(x) = . 

dx z 



4 This choice of the variables implies that t has the dimensions of a squared length and b of an 
inverse cube length, so that br 2 has the dimensions of a length. 
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Figure 2. Comparison between the Watson Airy function W(x) 
(red curve) and the ordinary Airy function (blue curve). 



The use of the unitary transformation 
(3.26) J(x) 
shows that it can be associated to the first ODE 



= ei a - W(x) 



(3.27) 



X 2 J(x) + l ( X ~ + ^ ) • /(.<") 



= . 



The transformation (3.26) has removed the second order derivative and has sim- 



plified the underlying group of symmetry, which has been reduced from SU(1, 1) 
to the dilatation group. We have mentioned this example because it may open 
a new interesting point of view to Schrodinger type equations involving quadratic 
potentials. 

Before concluding this paper we discuss whether the Airy polynomials can be 
exploited to obtain a series expansion of a given function. The problems associated 
with the orthogonal properties of ordinary and higher order Hermite polynomials 
has been thoroughly considered in a previous publications (see refs. [Tul lllj ): here 
we will use the point of view developed in ref. [11] . We consider indeed the following 
expansion 



(3.28) 



oo 

E 



a n Hi 3 \x,-\y\) , 



and we will show that the use of the operational tools developed in this paper pro- 
vides, in a fairly natural way, the coefficients a n , and we will display the orthogonal 
properties of this family of polynomials. On account of the definition of the higher 
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orders Hermite polynomials (see eq. (|l.6|)), eq. (3.28) yields: 



(3.29) 



The use of the identity ( 2.8 1 allows to specify the action of the exponential operator 
on the function f{x) as: 



(3.30) 



1 



d£ At 




E a " 



and the insertion of the explicit expression of the Airy function (see eq. (2.5)) 
yields: 



(3.31) 



1 1 



dt exp 



t 6 + i 



n=0 



The expansion of the a;-dependent part of the exponential on the lhs of the previous 
equation yields: 



(3.32) 



1 



2nn\ 3/3 



y\ 



dtf(t)d?Ai 



'3\V\ 



which holds only if the integral on the rhs of this equation converges. 

In these concluding remarks we have touched on several problems, which are 
worth to be explored thoroughly. The link between the Watson and the A^- 2p+1 ' > (x) 
functions, the possibility of obtaining more general Airy-type transforms and their 
relevance to Bell-type polynomials, and a more rigorous formulation of the or- 
thogonality properties of the Airy polynomials will be the topic of a forthcoming 
investigation. 
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